Matrix representation of the time operator 



Carl M. Bender°[j and Mariagiovanna Gianfreda'jj 

a Department of Physics, Kings College London, Strand, London WC2R 1LS, UK^ 
b Dipartimento di Fisica, Universitd del Salento and I. N.F.N. Sezione di Lecce, Via Arnesano, F73100 Lecce, Italy 

(Dated: January 19, 2012) 

In quantum mechanics the time operator Q satisfies the commutation relation [0,H] — i, and 
thus it may be thought of as being canonically conjugate to the Hamiltonian H. The time operator 
associated with a given Hamiltonian H is not unique because one can replace Q by + 0hom, where 
©hom satisfies the homogeneous condition [0hom,-H] = 0. To study this nonuniqueness the matrix 
elements of Q for the harmonic-oscillator Hamiltonian are calculated in the eigenstate basis. This 
calculation requires the summation of divergent series, and the summation is accomplished by using 
zeta-summation techniques. It is shown that by including appropriate homogeneous contributions, 
the matrix elements of simplify dramatically. However, it is still not clear whether there is an 
optimally simple representation of the time operator. 
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I. NONUNIQUENESS OF THE TIME 
OPERATOR FOR THE QUANTUM HARMONIC 
OSCILLATOR 

The time operator in quantum mechanics satisfies 
the commutation relation 

[Q,H] = i, (1) 

where H is the Hamiltonian [lj [2] . Even though time is 
not a dynamical variable, one may regard O as being a 
time operator because if one allows a normalized state 
\A) at * = to evolve to the state \B) = e~ im \A) at 
time t, then as a consequence of ([I]), 

(B\e\B)-(A\Q\A)=t. (2) 

It is not known how to solve (JlJ for the time operator 
for an arbitrary Hamiltonian H. Even for the compar- 
atively simple case of the quantum-harmonic-oscillator 
Hamiltonian 

H = \v 2 + \q 2 (3) 

the calculation of O is difficult because it involves singu- 
lar operators such as 1/p and 1/q. The time operator for 
the quantum harmonic oscillator was studied in detail at 
a formal level in Ref. [3]. The approach taken in [3] was 
to express O as a series of the form 

^ Q"m,n Tm , n ; (4) 
m,n 

where a m , n are numerical coefficients and T m ^ n (for inte- 
ger m and n) are a set of Hermitian basis operators. 
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For to, n > the operator T m n is defined as a totally 
symmetric average over all possible orderings of m factors 
of p and n factors of q. For example, 

To, o = 1, 

Ti,o = p, 

T i,i = \{pq + <iv), 

Ti,2 = l(pqq + qpq + qqp), 

T 2 ,2 = \{ppqq + qqpp + pqqp + qppq + qpqp + pqpq), 

and so on. The operator T mi „ is the quantum-mechanical 
generalization of the classical product p m q n . 

It is easy to evaluate commutators and anticommuta- 
tors of the operators T m rl . For example, as shown in 
Ref. [3], the operators T m n obey simple commutation 
and anticommutation relations: 

5 Tm,rt\ iTlT m n —i : 

[Q] -^m,n] — iTnTm— l,7ii 

\jPi -^"m,nj — ^T 7n ^-i n: 

{<7,??n,n} = 2T mj „_|_i. (5) 

By combining and iterating the results in ([5]), one can es- 
tablish additional useful commutation relations for T m n : 

[p : Tm,n\ — ^'^T m j r \ rL _\, 

[q 2 ,T m<n ] = 2imT m _ l!n+1 . (6) 

The totally symmetric operators T m n can be re- 
expressed in Weyl-ordered form [3]: 

Tm,n = ^±{ m X k q> m - k 

- (7) 

where m, n = 0, 1, 2, 3, • • • . The proof that T m ^ n equals 
the binomial-summation Weyl-ordered forms above is in- 
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ductive and requires repeated use of the Heisenberg alge- 
braic property that [q,p] — %■ The advantage of introduc- 
ing the Weyl-ordered form of T m ^ n is that it allows one 
to extend the totally symmetric operators T m ^ n either to 
negative values of n by using the first of these formulas 
or to negative values of m by using the second of these 
formulas. The commutation and anticommutation rela- 
tions in (5) and © remain valid when m is negative or 
when n is negative. 

In Ref. [3 it was shown that for the harmonic-oscillator 
Hamiltonian ^ the commutation relations in ^ may be 
used to obtain an exactly solvable recursion relation for 
the coefficients a m .„ in Q. The exact solution for has 
the series representation 



It is important to observe that the formal summation 



(-if 



fc+i 



^ 2k 

k=0 



'-2fc+l,-2fe-l- 



(8) 



In Ref. [3 the solution in ^ was called the minimal 
solution because of its relative simplicity. 

We emphasize here that m i n is merely a particular so- 
lution to the inhomogeneous linear equation ([I]) . To find 
other solutions to 0, we combine m i n with solutions 
to the associated homogeneous equation 



[h, e hom ] = o. 



(9) 



The subscript emphasizes that ©hom is a solution to the 
homogeneous equation ([9 1 . 

In Ref. [3] it was noted that the solution to ([!]) is not 
unique and it was claimed that any two operators satis- 
fying this commutation relation must differ by a function 
of the Hamiltonian H because this difference commutes 
with H. (We will see in Sec. [H] that this claim was actu- 
ally not correct.) 

One can interpret ^ by going to the classical limit 
— s- p m q n . In this limit the series ([sj) reduces to 

laical = £ ^ftr^+V 3 *- 1 = arctan^. (10) 
t 2k + 1 1 



Tn 



in (12 1 produces a bounded operator even though the 



individual terms in the summation are unbounded op- 
erators. (The same effect can be seen in the expan- 
sion e lq — J2^=o i n Q n /n\.) The fact that infinite sums 
of ill-behaved objects can be well behaved will be used 
throughout this paper. 

Observe that if the Hamiltonian H is symmetric under 
parity reflection [H, V] = 0, then as a consequence of 
|TJ the time operator O will also be symmetric under 
parity reflection. Also, if H is symmetric under time 
reversal [H, T] — 0, then O will be antisymmetric under 
T because the time-reversal operator changes the sign of 
i. (Note that if H is non-Hermitian but VT symmetric, 
then O will be VT antisymmetric [7].) The Hamiltonian 
for the quantum harmonic oscillator is both V and T 
symmetric because under parity reflection p — > —p and 
q — > — q and under time reversal p —> —p, q — > q, and 
i — > —i. Thus, the minimal solution in (Jsj) and (12) 



exhibits the expected V symmetry and T antisymmetry. 

The question addressed in this paper is similar to that 
addressed in Ref. [51 E] in which the uniqueness of the C 
operator in VT quantum theory was examined [7]. Be- 
cause the time operator is not unique, our objective here 
is to determine whether there is a solution to that is 
optimal in some sense. Our approach will be to calcu- 
late matrix elements of the time operator O and to ask 
whether there is a choice of homogeneous solution ©i lom 
which, when added to O m i n , gives particularly simple re- 
sults. 

In Sec. [n] we construct series representations for the 
solutions to the homogeneous commutator equation (Joj) . 
Then in Sec. |III| we show how to use zeta-function regula- 
tion to express matrix elements of the singular operators 
in these series representations as polynomials. Next, in 
Sec. IIVI we discuss the detailed mathematical structure 
of the polynomials obtained in Sec. |III| Finally, in Sec. [V] 
we give some concluding remarks. In the Appendix we 
explain the details of the zeta-function regulation tech- 
niques used in Sec. |III| 



Evidently, ©classical is conjugate to the harmonic- 
oscillator Hamiltonian in the sense of classical action- 
angle variables. 



One can also interpret 
using the simple identity [5 



8j) at the quantum level by 



1 (1 



2 \q 



(11) 



This formula allows us to sum the series in (|8| at the 
operator level (where the ordering of operators is crucial): 



\k+l 



mi „ = EyW(p 2fe+1 T 2/ 1 



4fc 

k=Q 
1 

-arctan I p 
2 Yq 



g-2*-yfc+l 



1 f 1 

-arctan —p 

2 \q F 



) 

(12) 



II. HOMOGENEOUS SOLUTIONS TO EQ. ([T] 
FOR THE HARMONIC OSCILLATOR 



In this section we show how to obtain homogeneous 
(nonminimal) solutions to the commutator equation ([T]) 
for the harmonic-oscillator Hamiltonian We will seek 
solutions of the general form Q. These solutions can be 
expressed as a single (not a double) sum on the index k 
and can be labeled by an integer parameter 7. 



A. Homogeneous solutions odd in p and odd in q 

In this paper we are only interested in homogeneous 
solutions that are odd in p and odd in q. The fact that 
H in (|3| is symmetric under space reflection implies that 
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cannot be odd in p and even in q, or even in p and 
odd in q. Furthermore, since H is symmetric under time 
reversal, must be either real and odd in p and odd 
in or else imaginary and even in p and even in q. In 
the latter case, however, will just turn out to be a 
function of H, and therefore the matrix elements of in 
the harmonic-oscillator-eigenfunction basis will be diag- 
onal. The point of this paper is to use the homogeneous 
solutions to simplify the matrix elements of the time op- 
erator, and therefore even-even solutions will not be of 
any value. Thus, we focus here on odd-odd homogeneous 
solutions. Solutions having this symmetry property can 
be written in the general form 



ft (7) _ n^T 



l,2fe+l- 



(13) 



k— — oo 



Substituting into ([£} gives the recursion relation 
(2fe + 1 - 2 7 )4 7) + (2fc + 3)4^ = (14) 



for —00 < k < 00. Thus, 



,(7) 
l k 



,/,"=(-l) fc oW 



0Fr(fc- 7 + l/2) 
2r(l/2 - 7)r(3/2 + k) ' 



(15) 



where is an arbitrary constant for each value of 7. 

The totally symmetric operators T m n that contribute 

to ©hom are shown in Fig. 
dots lie on a diagonal line. 
®min lie on a bold diagonal line. 



T 



as dots. For each 7 these 
he dots that contribute to 



To evaluate the sum in ( 13 ) for positive integers 7, we 

first use the symmetry property = a^fe-1+7 to split 
the sum into two sums. Then, we use the anticommutator 
relations in Q to obtain 



e 



(7) 

horn 



47 



fc=0 



^2 a k {?>••• U, U, T 2k+l ,-2fc-l}} • ■ ■} 



(27) times 

7-1 



+ {?,■■■ {P, {P, T-2k-l,2k+l}} ■ ■ •} 



(27) times 



a k J 2fc+l,-2fc-l- 



2 7 - 



(16) 



fc=0 



Finally, we use the identity ( 11 ) to simplify the totally symmetric operators of the form T_ m m . (This identity is valid 



for both positive- and negative-integer m.) This gives the result 

2fe+l 



o(7) _ 1 V fl (7) 
° l — — "27+1 



'horn 2 2 7 +1 ' J k 

fe=0 



2k+l 




^2 a k ^2fe+l,-2fe-l+27, (17) 



which we emphasize is only valid for positive-integer 7. 
It is possible to sum each of the four infinite series in ( 17 ) 
in terms of a hypergeometric function 



E 

k=0 



a 2* 1 



7,1, 



3 



(18) 



but there is no simple way to sum the finite series in ( 17 1 
The coefficients in 



15| are also valid for negative- 
integer 7, but now because 7 is negative, we use the com- 
mutator relations in ([5]) (instead of the anticommutator 
relations) to obtain the general formula 



m— 27,771 



(-!)>*-!)! 
(m + 27 - 1)! 
■ • [[T- mym ,q] ,q] ■ 



(27) times ' 



(19) 



and then we proceed as above. 

It is easy to see that Q^lL is not a function of H only. 

This is because the operator Q^om cnan g es sig n under 
time reversal T while H does not change sign. Further- 
more, if ©hom were a function of H only, then its ma- 
trix elements in the harmonic-oscillator-eigenstate basis 



would be diagonal. We will see in Sec. Ill that the matrix 
elements are not diagonal. 



B. Homogeneous solutions odd in p and even in q 

In this paper we are only interested in homogeneous 
solutions that are odd in both p and q, but we point out 
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FIG. 1: Schematic representation indicating the totally sym- 
metric operators T m , n that contribute to the sum in Q for 
m in and to the sum in |l3| for B^^. Each nonzero term in 
these sums is indicated by a dot. For each value of 7 these 
dots lie on diagonal lines. The dots associated with m i n lie 
on a bold diagonal line. 



that there are also solutions to ([9| having other symmetry 
properties. For example, solutions that are odd in p and 
even in q have the general form 



e 



(7) 



^ a fe 7 ''T , 2 7+ l-2fc,2/c, 



(20) 



where the parameter 7 = 0, ±1, ±2, .... 

Substituting ([20j) into we obtain the following two- 



(7). 



term recursion relation for the coefficients a 



4 7 +i(^+l)-(7-fc+l/2)ai 7) =0 (k = 0,1,2,...). (21) 



Unlike the recursion relation ( 14 ) , this recursion relation 



(7) 

is self-terminating; that is, if we choose a_{ = 0, then 
(7) _ n u ^ n Ji) is arbitrary, and a[7 ] ^ for 



a}." = for k < 0, a, 
(7) 

k > 0. In terms of the solution to this recursion 
relation is 



,(7) 



■ffn ^gr^.'ff <*=o,i,2,...). 



fc!r(-7- 1/2) 



(22) 

The series ( 20 ) with coefficients ( 22 ) can be summed 



as a binomial expansion: 



00 

k=0 



,(7) 



(1 



2N7+1/2 



(23) 



Thus, for 7 > the odd-even one-parameter family of 
solutions to the homogeneous equation Q is 




p p 



For example, for 7 = 



(27+I) times 



4 \ v p p v P p V V p p 



For 7 < the result is similar, but it involves commuta- 
tors. For 7 = — 1, we get 



9 



(-1) 

horn 



,(-!) r 



2/ 



p, arcsinh q 



arcsinh 



(26) 



Sec. 



while 6hom commutes with the Hamiltonian if in 



We can see again from (24) -(26) that, as we stated in 



([3]) , it is not a function of H . We verify this by observing 

that Qfo m is odd under parity reflection V while H is 
even. Moreover, in the classical limit, where p and q 



commute, (25) simplifies to 



«<°) 



e 



horn 



sgn(p)v / 2iJ. 



(27) 



The solutions in |24j)-(27l are solutions to the commu- 
tator equation ([9]), but they are physically unacceptable 
because they violate time-reversal symmetry. Thus, they 
must be rejected. However, the homogeneous solutions 
in ( 13 1, whose numerical coefficients are given in ( 15 ), are 



physically acceptable because they exhibit the correct be- 
havior under time and space reflection. This raises the 
question of which solution for the time operator is op- 
timal and in what sense it is optimal. We address this 
issue in the next section. 
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III. MATRIX ELEMENTS OF THE TIME 
OPERATOR 

We showed in Sec. |TT] that there are many possible 
choices for the time operator because there are an in- 
finite number of one-parameter families of homogeneous 
solutions to ([9| . One criterion for deciding on the optimal 
choice for the time operator (if there is an optimal choice) 
is mathematical simplicity. To investigate this criterion, 
in this section we calculate the matrix elements of the 
time operator in a basis consisting of harmonic-oscillator 
eigenfunctions. We begin by calculating the matrix ele- 
ments of the minimal solution (JsJ) . Then, we ask whether 
these matrix elements simplify when combined with ho- 



mogeneous solutions of the type in (13). 



In coordinate space the nth normalized cigcnfunction 
of the harmonic oscillator is 



ipn{<l) 



1 



7r 1 /4 A /2"n! 



e-« 2 / 2 H n (q). 



(28) 



If we try to calculate the matrix elements of Q min or 
@hom in this basis, we run into difficulties because each 
term in the series representation for these operators is 
singular (because powers of 1/p or l/q appear.) However, 
in this section we propose a method for circumventing 
these difficulties and thus for obtaining a finite result for 
the matrix elements. Specifically, we argue that while 
individual terms in the series representation for the time 
operator are singular, the sum of the series is a bounded 



operator [see, for example, (12)]. We then show that the 



(m, n) matrix element of a given singular operator in the 
sum can be made finite depending on the evenness or 
oddness of m and n and also depending on whether the 
singular operator is taken to act to the left or to the right. 



A. Matrix representation of the minimal solution 

To calculate the (m, n) matrix element of O m ; n , we 
must calculate the matrix element of each individual term 
^2fc+i,-2fc-i m the series |8]). The parity of this operator 
is even, so only the nonzero matrix elements correspond 
to values of m and n that are both odd or both even. 
Referring to (11), we can see that in coordinate space 
the (m,n) matrix element of ?2fe+i,-2fe-i has the form 



1 



H) 



2k+l 



2fe+l 



2k+l 



V>m(?) 9 q 



2fc+l 



Ml)- (29) 



Note that the first integral exists if n is even and the 
second integral exists if n is odd. When n is even, we 
reinterpret the differential operator in the second integral 
as operating to the left, and when n is odd we reinterpret 
the differential operator in the first integral as operating 
to the left. (Note that when the operator acts to the left, 
there is an additional minus sign because 2k + 1 is odd.) 



The only nonzero matrix elements of O m i n are 
(2n-2j|6 min |2n) 

1 ... fl 



y 1 

^ 2k + 1 



2n-2j| 



fc=0 

(2n-2j-l|6 min | 2ti-1) 



-a 



2fe+l 



E^n-^-ll^i) | 2 



\2n), 



2k+l 



(2n|6 min |2n - 2j) = -(2n - 2j|e min |2n), 
(2n-l|e roin |2n-2j-l) 

= -(2n-2j-l|6 min |2n-l), 



(30) 



where j — 1,2,3... and n > j. Note that the diago- 
nal matrix elements vanish: (m|0 m ; n |m) — 0. In Fig. [2] 
the nonzero matrix elements are shown as dots and the 
vanishing matrix elements are shown as zeros. 



The calculation of the matrix elements in ( 30 ) is 
lengthy, and we have relegated the discussion to the Ap- 
pendix. Here, we merely report the results of the calcu- 
lation. All the even-even and odd-odd matrix elements 
of the minimal solution (|8| can be expressed compactly 
as 

(2n-2j|6 min |2n) = l (-l) J 2-''/ 2 -i(j - l)! 

/ n!(2n-2j-~T)!i~ 
X y (n-j)!(2n-l)H W ' 

(2n + l-2j|e min |2n + l) = t(-iy +1 2^ 2 -\j - 1)! 

x./ ^-y + 1 »V o), (3D 



where 



(ra-j)!(2n + l)!! 



2P 



f(i) = S rw^i) ^' 



z(p) 



1 LV2+1J 



S P ,i v^U s. 



=1 p ] 



i - 1 

2r- 1 



o2r-2 "2r 



(32) 



(33) 



£?2r are Bernoulli numbers, and Sp i are Stirling numbers 
of the first kind. Table |T] gives the first eight values of the 
functions Z(p) and F(j). 

The explicit values of some matrix elements are then 



(2n-l|9 min |2n+l) 
(2n - 2\e min \2n) 

(2n-3|9 min |2n + l) 
(2n-4|6 min |2n) 



(2n + l)' 

2n 
(2n - 1) ' 




n{n — 1) 



(2n+l)(2n-l)' 

n(n — 1) 
(2n- l)(2n- 3)' 
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/OO. 0.0. 0.0.0. 
000.0.0.0.0.0 
.0 00.0.0.0.0. 
0.000.0.0.0.0 
.0.000.0.0.0. 
0.0.000.0.0.0 
.0.0.000.0.0. 
0.0.0.000.0.0 
.0.0.0.000.0. 
0.0.0.0.000.0 
.0.0.0.0.000. 
0.0.0.0.0.000 
.0.0.0.0.0.00 



Matrix representations of the homogeneous 
solutions 



The question is whether we can simplify the result 



in (31) by including the homogeneous solutions in (13 1 



} 



Let us first examine the contribution from the homoge- 
neous solution corresponding to 7 = 0. We find that the 

nonzero matrix elements for 6^ m lie on alternating diag- 
onals in Fig. [2j that is, the matrix elements are nonzero 
on the first dotted diagonals immediately above and im- 
mediately below the main diagonal. On the next dotted 
diagonals immediately above and below these dotted di- 
agonals, the matrix elements vanish; on the next dotted 
diagonals immediately above and below these dotted di- 
agonals, the matrix elements do not vanish; and so on. 
Specifically, we find that 



FIG. 2: Schematic view of the matrix elements of the minimal 
time operator m i n . The nonvanishing matrix elements are 
shown as dots. 



(n - 4j 



■2|e£Ll 



= -2<4 0) <n-4j + 2|6 r 



(34) 



3 


m 


P 


Zip) 


1 


1 


1 


1/2 


2 


2/3 


2 


-1/12 


3 


5/18 


3 


1/36 


4 


11/135 


4 


-1/90 


5 


101/5400 


5 


1/200 


6 


593/170100 


6 


-113/45360 


7 


2623/4762800 


7 


29/21168 


8 


383/5103000 


8 


-47/56700 



TABLE I: First eight values of the functions F(j) in (|32| and 
Z(p) in d33l. 



for j = 1,2,3,... and n > 4j - 2. 
which the matrix elements of 



(0) 



On the diagonals on 
do not vanish, each 
matrix element is a multiple of the matri x el ements of 



horn 



(0) 

'0 



= 1/2 in (151 and can- 
Thus, 



Hence, we can choose a 
eel half of the nonzero matrix elements of O m m 
including this homogeneous solution causes a dramatic 
simplification of the matrix elements of the time opera- 
tor! 



Next, we consider the matrix elements of O 



(7) 

horn 



for 



7 > 1. Here are some useful properties of the matrix 
representation of the operators T mj „: 



(2n-5|e mln |2n + l) 



hi / 2n(ra- l)(n- 2) 
9y (2ra + l)(2n- l){2n-Z)' 

(2n-6|e roin |2n) = 



hi / 2n(n-l)(n-2) 
~ 9y (2ra- l)(2n-3)(2n-5)' 
(2n-7|O min |2n + l) = 

44i / n{n - l)(n - 2){n - S) 



45 y (2n + l)(2n - l)(2n - 3)(2n - 5) : 
(2n-8|6 min |2n) = 

44z / n(n- l)(n - 2)(n - 3) 

45 V (2n-l)(2n-3)(2n-5)(2n-7)' 



(n - Aj\T. m ^ k \n) 
(n - 4j\T k k \ri) 



-{n - 4j\T k . m \n), 

0, (35) 



for n > 4j and for k,m > 1. For fixed 7, only the 
columns with j < [^r^J are not empty. This means that 



A (7) 



for 



the nonzero matrix elements of the operators Oj lom 
7=1 and 7 = 2 lie on just two diagonals above and two 
diagonals below the main diagonal; the nonzero matrix 
elements for 7 = 3 and 7 = 4 lie on two and six diagonals 
above and below the main diagonal; the nonzero matrix 
elements for 7 = 5 and 7 = 6 lie on two, six, and ten 
diagonals above and below the main diagonal; and so on. 
(The diagonals containing nonzero matrix elements are 
indicated in Fig. [3] by dots.) To summarize, 
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/ . 0000000000 
000.000000000 
.000.00000000 
0.000.0000000 
00.000.000000 
000.000.00000 
0000.000.0000 
00000.000.000 
000000.000.00 
0000000.000.0 
00000000.000. 
000000000.000 
0000000000 . 00 

V; ;;;;;;;;;;; ; 



/ . 000.000000 
000.000.00000 
.000.000.0000 
0.000.000.000 
00.000.000.00 
000.000.000.0 
.000.000.000. 
0.000.000.000 
00.000.000.00 
000.000.000.0 
0000.000.000. 
00000.000.000 
000000.000.00 

V : i 



I "horn ~ 



^0 0. 000. 000. 00 
000.000.000.0 
.000.000.000. 
0.000.000.000 
00.000.000.00 
000.000.000.0 
.000.000.000. 
0.000.000.000 
00.000.000.00 
000.000.000.0 
.000.000.000. 
0.000.000.000 
00.000.000.00 



FIG. 3: Schematic representation of the matrix elements for homogeneous solutions O^ om for 7 = 1, 2, . . . , 6. The dots indicate 
the nonzero matrix elements. For 7=1 and 7 = 2 there are only two diagonals where the matrix elements are nonzero. When 
7 = 3,4 there are four diagonals of nonzero matrix elements, when 7 = 5,6 there are six diagonals, and so on. 



(n-4j + 2|9 



(2j+o-l) 

horn 



2j+a-2 

= E «f' + ^ 1) («-4j+2|T 2fc+1 , 2(2j+CT _ fc _ 1) _ 1 |n) 



(36) 



fc=0 



where we have substituted 7 = 2j + a — 1 with <j = 0, 1, 2, The coefficients are given in 

The matrix elements of the operators T mj7l are associated with interesting classes of orthogonal polynomials. In 

particular, we are now interested in the case m, n odd. The jth class of orthogonal polynomials l'P^(n)} is 

I J k— 0,1... 

defined by the equation 



where 



(n - 2j\T 2k +i,2e+i\n) 



1 (n-2j)\ e ' k l+k - ]+ 
(-l) k j{£ + k)\{2k + 2)\{2£ + 2)\ 



i(n) (k,£>0,j>l,n>2j) 



2 e+k+2(i + + k + j + 1)\(£ + k-j + l)!(fc + 1)! ■ 



(37) 



(38) 



The polynomials "P^\n) (j = 1, 2, . . .) have degree k and argument n and are discussed in detail in the next section. 
Substituting (37) with j — > 2j — 1 and £ = 2j + a — k — 2 into (36), we obtain the following expression for the 

. — . . . . (-y) 

nonvanishing diagonals of elements in the matrix representation of the homogeneous solutions Q^Jm- 



(n-4:j + 2\e 



(2j-l+a) 
horn 



n) — —1a 



. ( a +23-i) (2j-l)(2a + 4j-2)! 



2 2i+<7-i( cr + 4 j _ 2)10-! V („_ 4j + 2)! 



(39) 



where j = 1, 2, 3, . . ., n > ij — 2, and a = 0,1, 2, 



IV. POLYNOMIAL REPRESENTATIONS OF 
THE HOMOGENEOUS SOLUTIONS 



where V^\x) — 1. The generating function 



In this s ectio n we describe the properties of the poly- 
nomials in (37). The polynomials V^\x) are orthogonal f or these polynomials is 
and satisfy the three-term recursion relation 



fe=0 



V^ +1 {x) = ^2x-2j + l)V^(x) 
+^(k + 2j)V^ 1 (x), 



>(J), 



nU)(t) - 02.7+1 ( 2 + l T 1 
° lj (2 



(41) 



(42) 



(40) 



From this generating function we identify these polyno- 
mials as Meixner polynomials of the first kind [51 [TP] . 



whose definition and properties are elucidated in the next 
subsection. In this paper j e N + , but if we set j = 0, 
we obtain the Hahn polynomials that were studied in 
Ref. QT]- 

The polynomials V k \x) do not have parity. However, 
if we shift the argument by x = iy + (2j — l)/2 and 

define Q ( k j) (y) = (-i) k V { k j} [iy + (2j - l)/2], then the new 
polynomials do have parity. Below we give the first six 
polynomials Q k \y) for j = 1 through j = 5: 



Qi 1] 



(.</) 



Q ( o 2) (y) 
Q { 2 \y) 



Q { o\y) 
Q ( 2 \y) 

Qf{y) 



= y 
= y 4 
= y° 



= i, 



Q ( o\y) 
Q?\y) = y 2 

Q [ :\y) = y 4 

y 5 



Qt\y) 



Q { i\y) 



3 
4' 

13 2 

2 y 
25 , 



= i, Q ( i\y) = y 



Q ( s\y) 

45 

16' 
309 

16"' 



5 

4' 

19 

— 1 
2 1 

35 



105 

16"' 
649 

16 



y, 



Qf\y) </ 
Q^(y) 



7 (V .„ 

4' 



25 



Q?Hv) = y 5 -'4y 3 + 



189 

16"' 
1109 



16 



Si 4) (2/) = V, 



9 

4' 
31 

7 

2 1 
55 



297 



16 ' 
1689 



16 



n 



Q? } (v) = y 3 -^y 



23 



29 



= 2/ 



(43) 



(44) 



(45) 



(46) 



for r 7^ s. We use the notation 

1$ = / dyw^(y)y 2 \ 



(49) 



to represent the moments of the weight functions w^\y). 
These moments are closely related to the Euler numbers 
E n . In particular, for the case j = (see Ref. [6]) 



(o) 



£ 2 fc|4- 



(50) 



In Table [TI| we list the values of the moments for j ranging 
from 1 through 6 and for k ranging from through 6. 



The moments defined in ( 49 1 and displayed in Table 
|IT| have some remarkable mathematical properties. To 
list these properties we first express these moments as 
polynomials IZk(j) = 4 fc ^^/(2j + 1) of degree k in the 
argument j. The first six of these polynomials are 

ftiC?) = i, 
n 2 {j) = 6j + 5, 

K 3 ( 3 ) = 60/ + 120j + 61, 

lZ A (j) = 840/ + 2940/ + 3486j + 1385, 

TZ 5 (j) = 15120/ + 80640/ + 163800/ 

+148800j + 50512, 
n 6 (J) = 332640/ + 2494800/ + 7595280/ 



+11641080/ + 8910726j + 2702765. 



(51) 



We then observe that these new polynomials resemble 
the Bernoulli polynomials and the Euler polynomials in 
that we recover the Bernoulli numbers B 2k and the Euler 
numbers E 2k when we evaluate the polynomials at special 
values of j. As we saw in (50), we obtain the Euler 
number \E 2k \ when we evaluate TZk(j) at j = 0: 



ftfc(0) = \E 2k \ 



(52) 



We obtain the Bernoulli numbers when we evaluate 
K k (j) at j = -1/2: 



H k (-l/2) 



2k y 



(53) 



Q { o\y) 
Q { 2 \y) 

Q { !\y) 

Qf\y) 



Qf\y) 



y- 



y 



n 

T' 

37 

~2 l 
65 



Qf{y) 

429 

"16"' 
2389 



-y 



16 



-y, 



35 



(47) 



The polynomials Q k \y) are orthogonal on (—00,00) 
with respect to an even, positive- weight function w^'(y): 



(48) 



There are many other values of j for which the polyno- 
mials TZk(j) give simple results: 

K k (-3/2) = (-4/- 1 , 
Kk(-l) = 

n k (i/2) = L^( 4 fc+i-i) 4 fe , 



K k (-5/2) = (-l)*- l 2 



fc-l 2fc-3 /r>2fc— 2 



n k (-2) = 



(54) 



9 



k 


4Vi 1} 


4Vi 2) 


4Vi 3) 


4^i 4) 


4Vi 5) 


4Vi 6) 





1 


1 


1 


1 


1 


1 


1 


3 


5 


7 


9 


11 


13 


2 


33 


85 


161 


261 


385 


533 


3 


723 


2705 


6727 


13509 


23771 


38233 


4 


25953 


134185 


426881 


1045161 


2170465 


4024553 


5 


1376643 


9451805 


37611847 


110683809 


268614731 


570100453 


6 


101031873 


892060285 


4355312801 


15209937261 


42750877345 


103164046973 



TAB LE II: Numerical values of the moments /ij. defined in 
(49l. 



A. Connection between Vjf (x) and 
Meixner /Meixner-Pollaczek polynomials. 



The polynomials Vjf' (x) defined in (37) satisfy the 



three-term recursion ( 40 ) and are related to the Meixner 
polynomials: 



vi j \x)=2- k M k (x-2j, 2J + 1.-1), 



(55) 



where the Meixner polynomials M n (x, b, c) belong to the 
classical orthogonal polynomials of a discrete variable 
and may be defined in terms of hypergeometric functions: 



M n (x, b, c) = T{ ^ k \ F 1 (-k, -x- b: 1 - c- 1 ) 
T(b) 



(56) 



For real values of the parameters b > and < c < 1 
they satisfy a discrete orthogonality relation with weight 
function 



w(x) 

whose moments are 



~(b,c) 



E 

n=0 



T(b + x)c x 
T(b)x\ ' 

T(b + n) c n n k 
T(b) n\ ' 



(57) 



(58) 



The problem with evaluating the moments associated 
with the polynomials V k \x), that is, the moments of 
the Meixner polynomials with c = — 1, can be solved by 
performing a zeta- function regularization of the divergent 
series in ( 58 1 : 



-(6-1) 



= lim > 



(-l) n T(b + n)n 
T{b)n\ 



k+s 



(M0). (59) 



The polynomials (%) are defined as 



Q<j\ x ) = (-i)« V U'[ix+(2j-l)/2}, 



(60) 



which are listed in (43 )-(|47j), satisfy the three-term recur- 
sion relation Q.{ j l{x) = x Q { k J \x) - |(fc + 2j) Q^(x), 



and are connected with the Meixner-Pollaczek polynomi- 
alsP* (*;</>) H2| 

Q k j) (x) = 2- k k\P [ k J+1/2 \x; it/2). (61) 
The Meixner-Pollaczek polynomials are defined as 

Pk(x > 0) = T{2 ^\)kT 2Fi ix + A; 2A; 1 - e_2i0) ■ 

(62) 

For < <f) < ir and A > 0, Pk(x, 4>) are continuous orthog- 
onal polynomials on the real line with respect to a contin- 
uous positive real measure given by the weight function 

w{x) = ^{2sin(j>) 2X \T{\ + ix)\ 2 e x{2 ^\ (63) 

The Meixner-Pollaczek polynomials P k (x) are analytic 
continuations of the Meixner polynomials M. n (x, b, c) in 
the parameter c |13j . The connection is given by 

p£( x ) = e - ik *M k (ix - A, 2A, e'^/kl, (64) 



which is like the relation between Q k \x) and V k '{x) 



>U)t 



B. Evaluation of moments associated with Q.\?\x) 



Equations (61 ) and (64) show the relationship between 



Q k 3 \x) andX fe (x): 

QkH*) = H) k 2- k M k (ix-j - 1/2, 2j + 1, -1). (65) 

Moments of the polynomials Q k (x) can be evaluated in 
two ways. One way is to use ( 64 ) to express the moments 

fjP^ associated with the polynomials Q k \x) in terms 



of the regularized moments in ( 59 ) associated with the 
Meixner polynomials 



(:;) 



H) fe E 



a=0 



0' + l/2) a 



xV™^^. (66) 
^ n!r(2j + l 



10 



For (2j + 1) > we can express the gamma function 
in (66) as a sum of powers of n, 



r(2j + i 



2j P /R 



'2j 

(7 = 



(67) 



and then the zeta regularization ( 59 ) gives 



o2j + l 2/c-l /oi, 

(-!)*{ (j + l/2f k + - ^ ,J ' 



r(2j 



2,; 



0- + l/2) a 



EC" 1 )^^ E J (2 2fe+ — +1 - 1) C(-2fc + « - a) } . 



(68) 



r 



Equation (68) reproduces the particular values of 

TZk(j) — V(2j + 1) given in this paper for j = 
and j = 1/2. The special case lZk(— 1/2) is obtained by 
first dividing ( 68 ) by (2j + 2) (2j + 1) and then substitut- 
ing the value j = —1/2 in the result ( [53|. 

For ( 2 J + 1 ) < the series over n in (|66[) becomes finite 
and regularization is not necessary. The result is 



(-l)^ 1 ^^) (i + l/2) a 

xr(-2i) ^ — — — — 



(69) 



which gives the same values for 1Z n (— 1), 7£ n (— 5/2) de- 
rived earlier in this paper. 

The other way is to evaluate the moments of the poly- 
nomials Q^\x) using their relation with the Meixner- 



Pollaczek polynomials (61). Using the result [14 



dxe- i7I - 2<t,)x \r(tx + a)\ 2 = 7rr(2a) (2 sin $y 2a , 

> 

(70) 

differentiating the above formula n times, and dividing 
both sides by 2™, we have 



dxx n e- (7T - 



2<h)x 



T{ix + a)\' 



2-Nrr(2a) — (2sin0)- 2a . 



(71) 



From (|7T|) with <j> = it/2 and a = j + 1/2 and using (61 ) 



and (|63|) we obtain the moments of Qu (x): 

,-2n-l u ^o„:„ m-2j-1 



(2 sine 



i=7r/2 



(72) 



This result is equivalent to ( 68 ) and ( 69 ) 



Favard's theorem [T5] states that if a sequence of monic 
polynomials satisfies a three-term recursion relation of 



the form xP n (x) — P n+ i(x) + a n P n (x) + b n P n -i(x) with 
a„,!)„ G K and b n > 0, then there exists a positive 
measure respect to which the polynomials are orthog- 
onal on R. In our case, because the recursion relation 

for the polynomials Q ( i\x) is xQ$(x) = Q^l^x) + 

n(n + 2j)Q^}_ 1 {x), they are orthogonal with respect to 
a positive measure for j > —1/2, in agreement with the 
results in (52) - (54). In Ref. [TB] an extension of the 



Meixner-Pollaczek polynomials for A < is studied and 
a nonstandard inner product is defined with respect to 
which they are orthogonal. 



V. CONCLUSIONS 

In this paper we have investigated the question of 
whether the minimal solution for the time operator can 
be significantly simplified by adding homogeneous solu- 
tions. To answer this question we had to make sense of 
and evaluate infinite sums of singular operators, and to 
do this we calculated matrix elements of these sums and 
then used zeta summation to sum the series. The ma- 
trix elements give rise to remarkable Meixner polynomi- 
als having interesting mathematical properties in which 
Euler and Bernoulli numbers repeatedly appear. 

Our principal conclusion is that the matrix elements 
of the minimal solution (n — 2j|O min |n) for j = 1,2... 
and n > 2j cannot all be canceled by adding one or more 
homogeneous solutions ©u™ for 7 = 0,1,2,... to the 
minimal solution. Specifically, for the minimal solution 
the matrix elements (n— 4j+2|0 m ; n |n) for j = 1, 2 . . . and 
n > 4j — 2 can be canceled by adding to O m i n the homo- 
geneous solution 0jjo* m with the choice of the parameter 



,(0) 



1/2. However, the homogeneous solutions 8 



horn 



for 7 > 1 cannot cancel all the matrix elements of the 
minimal solution (n — Aj + 2\<d m - m \n) for j = 1,2... and 
n>4j- 2. 

To elaborate, we have shown that for a finite value of 



7, the number of nonzero columns of the operator 8 



«(7) 



horn 



11 



is less than the number of the nonzero columns of 9 m i n . 
This means that it is impossible to cancel all the matrix 
elements (n — 4j + 2|0 min |rt). As 7 — > 00, the number of 

nonzero elements in the matrix representation of ©ulL 
becomes equal to the number of nonzero elements in the 
matrix representation of m m- This means that we can 
then consider the operator 



7=0 



(73) 



The matrix elements in the first nonzero columns above 
and below the principal diagonals (J — 1) are 



In - 2\Q\n) 



V n(n - 1) X] °o 



(g+ D (2a + 2)! 



(cr + 2)!2 CT er! 



If we choose 



( CT ) ((7 + 1)! ! 

l ° " (2a)! {2t) ' 



(75) 



with t to be determined, the matrix elements become 

(* — 1)* r~n 



(n-2\e\n) 



2t 



n — 1 



+4i v ^l) |±|^ , (76) 



where we have used the generating f unc tion of the poly- 
nomials Va (n) to sum the series in ( 74 1. The values of t 



that make the matrix elements ( 76 1 vanish are the roots 
of the equation 



t(t - 2)" +1 - 8(n - l)t - (t - 2)™ +1 =0 (n > 2), (77) 
with t + 1 0, 2. 

This demonstrates that to cancel all the elements in 
a fixed diagonal in the matrix representation (n — Aj + 

2|0|n), the choice of the parameter must be different 
for each n. In general, the value of that cancels the 
matrix elements (n — 4j + 2|0|n) depends on both j and 
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Appendix A: Calculation of Matrix Elements 

In this Appendix we show how to calculate the matrix 



elements in (31). Let us consider the operators 



l__d 
q dq 



Jk 



d 1 
dq q 



(A: = 1,2,...). (Al) 



Here, Ik acts on even eigenstates according to 



k k—a / , \ 

W = EE4« ^„)R nia , fi q- 2 P- a H n - a (q)e-^-H n (q)e-^ 3 , (A2) 



P 

where H n are the Hermite polynomials. Similarly, acts on odd eigenstates and gives the result 

k k—a 

y a + j3 



JkMQ) = E E ( 1 R ) Rn, a ,f>q~ W ~ a Hn- a e- x2/2 . (A3) 



a=0 /3=0 



The numerical coefficients R n , a ,f3 are given by The Hermite polynomials have the explicit form 



L(n-«)/2J 

(-l)°+i2°-g(q + 2fl!n! (n-a)l , 

Rn ' a>P ~ ^V^\a\p\(n-a)\ ' (A4) C7«,„,, - (-1) 2 z!( „_ a _ 2Z) ,- ( A5 ) 



12 



Thus, we can evaluate the integral 

/>oo 

S n , m ,u = / dqe-" 2 q n - 2 "H m (q) 



{ 2 m - 1 (n/2-v)\r(n/2-v+l/2) r < _ „ s 
_ J (n/2-i/-m/2)l - " LV >' 

I otherwise. 



where a mi 
(n — m)/2 — a — /?, and 



(n - m)/2, /3 max = (n - m)/2 - a, 4 



The matrix elements (m|Jfc|n) and (m|«7fc|n) can then 
be written in the form 



Omax /3max f m « 

>i« = E E E 

Q = l 0=0 £=0 
C*max /9max £ ma x 



a (2k + 1 

a + 2/3 Va + /3 



(2fc + 1 N 
„ ) U n , m , a ,t,p> (A6) 
■ a ^ ' 



7 r 1 /4^2 m m! 

(A7) 

The coefficients £/ have the form 



2n+l,2m+l,a,£,/3 



?7 : 



2n,2m,a,t,j) 



( — l) Q+£+1 2 2a+,9+m_n 



'(2n + l)! (a + 2/3)!(2n-2a-2^-2f + 1)! 

(2m + 1)! /3!#a!(n - m - a - f3 - £)l(2n -a-2l+l)V 



' (2n)! (a + 2/3)!(2ra - 2a - 2(3 - 21)! 

(2m)! /3!£!a!(n - m - a - ,8 - £)!(2n - a - 2f)! ' 

I 



With the change of variables n — 2n' and m = 2n' — 2j 
for j = 1, 2, ...n', we make the replacement 



U, 



n,j,a,e.,0 — U2ri,2n'-2j : a,£ 



(A8) 



and (A6) becomes 
(m|J fc |n) = E E E r fl )0».W,/J. ( A9 ) 

a=0 /3=0 f=0 ^ 1 ' 

Then after the second change of variables, a = p — a, 
P = a, U n ,i,p,l,a = U n ,j.p-a,iM, we obtain 

(m|/*|n) = E E E?V^( J^wA^ 
j p—l i-p 



HJ k \n) = E E E 

p=l cr = 1 = 

Next, we use the formulas 
(2n-2j|e min |2n) = 



2fe + l 
P 



(A10) 



2 E^XT< 2n - 2 J'lW|2n>, (All) 

fc=0 



(2n + l-2j|O min |2n+l) 



X)^q^<2rc+l-2j|W|2n+l), (A12) 



fc=0 



to rewrite the even-even and odd-odd matrix elements 
(2n-2j|8 min |2n) = 



2'' 



oo 1 

y 1 

f-^ 2fe + l 



(2n + l-2j|6 mln |2n + l) 



^° y 



(S-l)\(j-8)\ 



oo 

y 1 

^ 2fc + 1 



2fc + 1 
<5 



2fe + l 



(A13) 



where for n = 1, 2, . . . and j — 1,2, ... n 

' n\(2n- 2j - 1)!! 



ft! 



l -(-l)i 2 ^-i( i _i) ! 



{n-j)\{2n-l)\V 

1 n\{2n- 2j + 1)!! 
(ra- j)K 2 "+ 1 ) !! ' 



Now an apparent problem arises: The series over k in 
(A13) is divergent. To overcome this problem we resort 



to zcta summation to evaluate the divergent series. Zeta 
summation is conventionally used to evaluate a divergent 
sum over modes. (For example, it is used in the calcu- 
lation of the Casimir force [5].) The idea behind zeta 
summation is to extend £(s) = 53&=i ^ -S i which is valid 



13 



for Res > 1, to negative values of s. To do so, one must 
use the functional equation for the Riemann zeta function 

C(l -s) = 2(2^)-*r( S ) cos (^) C(«). (A14) 

We can thus express the divergent series YlkLi a (^) as 
lini s _ > o^^_ 1 a(k)k~ s . This gives results such as 

oo oo 1 

c(o) = ^fc° = -i c(-i) = £ fc = -7^ 



k=l 

In our case the sum 



fc=i 



12 



Z M "2^2fcTT 

fc=0 



becomes 



z ip) = S P ,i + limy" 

K — >- * ■ 



2fc + l\ k 



s-to ^— ' \ p / 2fc + 1 



(A15) 



(A16) 



It is now straightforward to evaluate the sum over k for 
generic values of the parameter p. We first rewrite the 
sum by expressing the binomial coefficient as a polyno- 
mial in the variable k 



2k + 1 
P 



p gi 



= J2^(2k + iy (p= 1,2,3...), (A17) 



where i are the Stirling numbers of first kind. Then. 



3(p) = *p.i+E%E( V) - (AIN) 

i=l r=0 



Because ((0) = -1/2 , ((- 2n) = and f(l - 2n) = 
—B2n/(2n) for n = 1,2,3..., (A16| is nonzero only for 
odd values of r and for r = 0. Replacing r — > 2r — 1, we 
obtain the final results in ( 32 ) and ( 33 ) . 



[1] E. A. Galapon, Proc. R. Soc. Lond. A 458, 451 (2002). 
[2] E. A. Galapon, R. F. Caballar, and R. Bahague, Phys. 

Rev. Lett. 93, 180406 (2004) and Phys. Rev. A 72, 

062107 (2005). 

[3] C. M. Bender and G. V. Dunne, Phys. Rev. D 40, 2739 
(1989). 

[4] C. M. Bender and G. V. Dunne, Phys. Rev. D 40, 3504 
(1989). 

[5] C. M. Bender and S. P. Klevansky, Phys. Lett. A 373, 
2670 (2009). 

[6] Another simple relation is T„ t „ = P n (Ti,i), where P„ is 
the nth Hahn polynomial; this was derived in C. M. Ben- 
der, L. R. Mead, and S. S. Pinsky, J. Math. Phys. 28, 
509 (1987). 

[7] C. M. Bender, Rept. Prog. Phys. 70, 947-1018 (2007). 
[8] K. A. Milton, J. Phys. A: Math. Gen. 37, R209 (2004). 
[9] A. Erdelyi, W. Magnus, F. Oberhettinger, and F. G. Tri- 
comi, Higher Transcendental Functions (Krieger, New 



York, 1981), Vol. II, pp. 224-225, 1981. 
[10] E. W. Weisstein, "Meixner Polynomial of the First 

Kind." From MathWorld-A Wolfram Web Resource. 

http://mathworld.wolfram.com 

/MeixnerPolynomialoftheFirstKind.html 
[11] C. M. Bender and G. V. Dunne, J. Math. Phys. xx, 3504 

(1988). 

[12] R. Koekoek, R. F. Swarttouw, The Askey- scheme of hy- 
pergeometric orthogonal polynomials and its q-analogue, 
Delft University of Technology and Systems, Department 
of Technical Mathematics and Informatics, Report No. 
98-17, (1998). 

[13] N. M. Atakishiyev and S. K. Suslov, J. Phys. A: 

Math. Gen. 18, 1583 (1985). 
[14] F. Pollaczek, C. R. Acad. Sci. 230, 1563 (1950). 
[15] J. Favard, C. R. Acad. Sci. Paris 200, 2052 (1935). 
[16] T. K. Araaya, J. Math. Anal. Appl. 305, 411 (2005). 



